* Classical statitistical mechanics (some resonable combination of ...)
** equiprobable microstates, and definition of entropy.
** Boltzman factor: P(E) proportional to e^(-E/kT).
** Equipartition theorem.
** Maxwell's velocity (and corresponding speed/KE distributions).
** P(E) = e^(-E/kT) / Z, Z = sum e^(-E/kT).

* Quantum statistics
** Reminder of symmetry requirements for bosons and fermions
** Corresponding changes in distribution functions
f(E) = 1/[Ae^(-E/kT) +/- 1] = 1/[e^(-(E-E0)/kT) +/- 1]
E0 = fermi level for fermions, E0 = mu (chemical potential?) for bosons
(from Liboff).  Physical meaning.

* Applications
** Blackbody radiation.
** Einstein's model for heat capacity (atoms all oscillate with equal
frequency; see Callan ch 15 for alternative analysis and ch 16 for
the Debye model).
** Liquid Helium (bose condensation) -- qualitative only.
** Free electron treatment of conduction within a metal.

I need to research this stuff.


Possible explanation of occupation factors (based on grand canonical
formalism):

mu = energy that a particle will have if it decides to be "elsewhere"

for fermions, two choices: either the particle is "elsewhere" with
energy mu or else it is in some particular state we are focussing on,
with energy E.  The relative probabilities of these two possibilities
is given by the boltzmann factor:
P_occ = e^(-(E-mu)/kT) P_unocc

put another way: P_occ = e^(-(E-mu)/kT) / Z and P_unocc = 1/Z, where
Z is a normalization factor: Z = 1 + e^(-(E-mu)/kT).  This reduces to
f(E) = P_occ = e^(-(E-mu)/kT) / [1 + e^(-(E-mu)/kT)] =
     = 1 / [e^((E-mu)/kT) + 1]

for bosons, there is no limit in how many particles can be in that
state.  assuming that bosons can have energy mu when "elsewhere"
implies relative probabilities (g = (E-mu)/kT)
e^0, e^(-g), e^(-2g), e^(-3g), etc when N = 0, 1, 2, ....  Absolute
probabilities can be obtained with help of normalization
Z = 1 + e^(-g) + e^(-2g) + ... = 1 / [1 - e^(-g)].
The probability for occupation number N is then
P(N) = e^(-Ng)/Z.
Now,
N_av = sum_N N P(N) = sum_N (N+1) P(N+1) = e^(-g) sum_N (N+1)P(N)
     = e^(-g) (N_av + 1)
Solving for N_av yields
f(E) = N_av = e^(-g) / [1 - e^(-g)] = 1 / [e^g - 1].

for photons, "elsewhere" means "does not exist", so mu=0, and we
get the blackbody spectrum (after considering density of states).

Note: if particles were to be treated as distinguishable, then
we would have combinatorial factors going along with each e^(-Ng).
That would have consequences.  I need more research.


>> A macrostate is a collection of microstates consistent with certain
bulk measurements (e.g., take a bucket of water: an example of a
macrostate would be specification of total energy, volume and number
of particles in the bucket, or perhaps some specific DISTRIBUTION of
energy throughout the volume of the water / perhaps E1 and E2
separately for sub-regions 1 and 2).  Entropy for a given macro-state
= k ln(# microstates consistent with that macrostate).

Fundamental principle of thermodynamics: At equilibrium, the observed
macrostate is the most likely one (i.e., highest entropy).  For large
systems, the probabilities are so skewed (factors of e^(10^23) are not
unexpected) that the lesser likely macrostates are NEVER observed.
Determining entropy as a function of internal energy, volume, total
number (etc) completely determines the equilibrium behavior of a
substance.  All thermodynamic quantities related to equilibrium can be
defined in terms of this function.

Proposed sequence of discussion:

Insert somewhere in this discussion: (1) equipartition, (2) maxwell
velocity (speed/KE) distribution.

* isolated system (fixed system energy E): all states equiprobable.
Entropy discussion above.

* system connected to thermal bath (E can vary): probability of
state i is e^(-E_i/kT) / Z, where Z = sum_j e^(-E_j/kT) is the
partition function (normalization).

P_i = e^(-E_i/kT) / Z; Z = sum_j e^(-E_j/kT)  [state focussed]
P(E)dE = g(E)dE e^(-E/kT) / Z; Z = integ g(E)dE e^(-E/kT)  [energy focussed]

generally, the most probable energy is NOT the lowest possible value,
since g(E) is usually very small there (the ground state as a single
state is the most likely state of any other state, but the ground
state ENERGY is unlikely to be the most likely ENERGY, because there
are so few states with that energy).

* if system consists of completely non-interacting particles, each
particle can be treated as its own isolated system: all particle
states (with fixed E) are equiprobable, but there are no transitions
(because no interactions), so in fact, the states remain fixed.  This
is completely unrealistic for large systems.

* if system consists of weakly (but not completely non-) interacting
particles (much more realistic assumption), then we may approximate
(at any one time!)  psi_system = psi_1 psi_2 ... and E_system = E_1 +
E_2 + ... (where 1, 2, ... are individual particles).  in this case,
particles make up their own system connected to a thermal bath.
probability that particle i is in state j is given by e^(-E_ij/kT) /
z_i, where z_i = sum_j e^(-E_ij/kT).

* if we focus on two particles (1 and 2) of our larger system, then
probability that particle 1 is in state i and particle 2 is in state
j is [e^(-E_1i/kT)/z_1]*[e^(-E_2j/kT)/z_2] = e^(-(E_1i + E_2j)/kT) / z_12
where z_12 = z_1 z_2.  Note that these probabilities can be read directly
from the combined partition function z_1 z_2 expanded term by term.
more than two particles handled similarly.  For weakly interacting
particles, the partition function factors: Z = z_1 z_2 ... z_N.

* suppose we have N particles, each can exist in (the same) two
states.  we are not so concerned with which particles occupy which
states, but only on how many particles occupy state A and how many
occupy state B.  Then z_N = z^N = (e^(-E_A/kT) + e^(-E_B/kT))^N
= sum_(nA,nB; nA+nB=N) (N!/nA!nB!) e^(-(nA*E_A+nB*E_B)/kT) (binomial
expansion).  It is clear that P(nA,nB) (probability that occupation
numbers are specifically nA and nB, nA+nB=N) is one of these terms
divided by the total partition function.  If we focus on just one
of these states
P(nA) = (N!/nA!(N-nA)!) e^(-(nA*E_A+(N-nA)*E_B)/kT) / z_N
      = (approx) C(nA) e^(-nA*(E_A-mu)/kT) / z_A
where C(nA) is a combinatorial factor, and mu is some sort of
"average" energy per particle that decides NOT to be in state A
(so E_A - mu represents the energy penalty associated with the
particle ending up in state A rather than elsewhere), and z_A is
the normalization factor appropriate for this form of the equation.
This approximation becomes more accurate as the total number of
particles (N) grows large, because P(nA) becomes more and more peaked
around a specific narrow range of values (so mu can be treated as
independent of nA).  For more than two states, P(nA) would require a
summation over the other occupation numbers, making mu and C(nA) more
complicated (actually, I believe C(nA) = N!/nA!(N-nA)! after all the
dust settles), but still a good approximation.

CHECK THIS: <nA>/<nB> = e^(-E_A/kT) / e^(-E_B/kT) (using exact P),
even with more than two states in the picture.  Of course, this will
go out the window when we consider bose/fermi statistics because the
combinatorial factors change.

* Above analysis is consistent with classical physics, and also with
QM if the particles are distinguishable (not identical).  Proper
quantum treatment of N identical particles requires some care.

In QM, the wave-functions for a system of identical particles must
obey proper symmetry (completely symmetric for bosons, completely
antisymmetric for fermions).  for two particles, one in state A and
one in state B, the 2-particle wavefunction looks like this:
psi(1,2) = psi_A(1)psi_B(2) +/- psi_B(1)psi_A(2) [normalize]
where + is for bosons and - is for fermions.  This can be generalized
for N particles.

This results in a change in the combinatorial factor C(nA).  In a
two-state system, rather than having N!/nA!nB! separate states
corresponding to state occupation numbers (nA,nB), there is at most
one state (the completely (anti-)symmetrized combination N-particle
state), and, if the particles are fermions and either nA or nB is > 1,
zero states.

If we focus on a single (orbital) state, A, then
P(A) = C(nA) e^(-nA*(E_A-mu)/kT) / z_A
where C(nA)=1 for bosons (all nA) and fermions if nA <= 1, and C(nA)=0
for fermions with nA >= 2.  The normalization, z_A, is given by
z_A(boson) = e^0 + e^(-(E_A-mu)/kT) + e^(-2(E_A-mu)/kT) + ... 
  =  1/(1-e^(-(E_A-mu)/kT)), or
z_A(fermion) = e^0 + e^(-(E_A-mu)/kT).

After a bit of math, the average occupation number, <nA>, is given by
<nA>(boson) = sum_nA nA P(nA) = 1/(e^((E_A-mu)/kT) - 1), or
<nA>(fermion) = sum_nA nA P(nA) = P(1) = 1/(e^((E_A-mu)/kT) + 1).
and the average energy of all particles in this state is
<E> = E_A <nA> = E_A / (e^((E_A-mu)/kT) +/- 1).

* Note that these average occupation numbers
<nA> = f(E_A) = 1/(e^((E_A-mu)/kT) +/- 1)
is for one given state.  The number of particles within a given energy
range is given by
N(E)dE = D(E)dE f(E),
where D(E) is the density of states, which are approximated as a
continuum (normally a reasonable approximation for large systems).

* Bose condensation: treatment of states as a continuum puts a strict
upper bound on the quantity N lambda^3 / (2S+1)V of 2.612 (=sum
n^(-3/2)) (lambda = hbar / sqrt(2pi m kT) = "thermal wavelength") [ch
18 Callan].  Big deal, you say?  Well, when T becomes too small, this
bound is violated.  The reason for the violation is that as T->0,
mu->E_gnd from below, and causes f(E_gnd) to diverge.  In this case,
treating the ground state as part of the continuum breaks down because
its occupation number grows out of proportion to the fact that it is
only one state out of (on the order of) 10^23.  Proper treatment of
this state reveals the large occupation number, which continues to
grow as the temperature continues to fall: <n_gnd>/N = (1 -
T/Tc)^(3/2), where Tc is the critical temperature where Bose
condensation occurs.  This phase transition is also marked by a sudden
jump in the specific heat [p. 420 of Callan has some graphs]

The particles occupying the ground state behave in a totally
non-classical way: the motions are totally corrollated, there is very
little dissipation possible (those particles can't lose energy: they
are in the GROUND STATE), so they essentially move without friction
(non-viscous flow).  This accounts for the superfluidity of He_4, and
also serves as an explanation of superconductivity (electron (Cooper)
pairs form at low T, and they act as bosons), as well as the eventual
superfluidity of He_3 (a fermion, but, like electrons in a
superconductor, they too can form Cooper pairs at extremely low T).
